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Abstract: In this paper analytical results are presented for higher order corrections to 

coefficient functions of the operator product expansion (OPE) for the correlator of two 

pseudoscalar gluonium operators 0\ = G^G^. The Wilson coefficient in front of the 

scalar gluon condensate operator 0\ = —\G^ U G^ V is given at three- loop accuracy. The 

leading coefficient Co in front of the unity operator Oq = 1 has been calculated up to 

l^j three-loop order some time ago [1] but has been checked independently in this work. It is 

,/ ■ interesting to see that the coefficient C\ in the pseudoscalar case is finite, whereas contact 

CD ■ terms appear in Cq in this case and in both coefficients Cq and C\ in the cases of the scalar 

gluonium correlator and the energy momentum tensor correlator [21. For the corresponding 

Renormalization Group invariant Wilson coefficients which are also constructed the results 

are partially extended to four-loop accuracy. All results are given in the MS-scheme at 

zero temperature. 
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1 Motivation 

Euclidian correlators of local operators are important objects in quantum field theory. 
Firstly, they have many important applications, e.g. in sum rules, where they are connected 
to physical quantities like spectral densities through dispersion relations. Secondly, they 
often have interesting properties in themselves, like their non-trivial renormalization, which 
are important for the understanding of quantum field theories. Such correlators are defined 
in momentum space as 

i [d*xe iqx T{[O}(x)[O}(0)} (1.1) 

with a large Euclidian momentum q. Here and in the following the squared brackets indicate 
that the renormalized form of some operator O is used. Usually, we are interested in the 
vacuum expectation value (VEV) of the correlator 

n(Q 2 ) = i [d A xe^ (0|T{ [O](x)[O](0)]}\0) (Q 2 = -q 2 ) (1.2) 



which can be calculated in perturbation theory. But if we take |0) to be the physical vacuum 
state we also have to consider non-perturbative effects. Starting from the perturbative 
region of momentum space this is done by means of an operator product expansion (OPE). 
The idea is to expand the bilocal operator product (1.1) in a series of local operators with 



Wilson coefficients depending on the large Euclidean momentum q [3]: 1 

/ . ,^ „ 2 dim(0)-dim(O i )-4 „ 
d 4 xe^r{[o](x)[0](o)} = £cf( g )(Q 2 ) a qB (L3) 

i 

= £G^)(Q 2 ) a [OJ- (1-4) 

i 

2 dim(Q)-dim(O i )-4 

where the index B marks bare quantities and the factor (Q ) 2 constructed 

from the mass dimensions of the operators involved makes the Wilson coefficients Ci{q) 
dimensionless. 

In a sum rule approach to glueballs three operators are usually investigated as insertions 
on the lhs of (1.3) (see e.g. [5]): 



Ox(x) = —GTG^x) 




(scalar) 


(1.5) 


6 1 (x) = G^G^{x) 




(pseudoscalar) 


(1.6) 


0^{x) = T» v (x) 




(tensor) 


(1.7) 


where G^ u is the gluon field strength tensor, 








G^y = £ 


G pa 




(1.8) 



the dual gluon field strength tensor and T^ u the energy-momentum tensor of QCD. Having 
discussed the correlators of Q\ and Oj? in [2] the results for the correlator of (1.6) 



X t (q) :=tydW«*r{[6i](x)[6i](0)}, (1. 



whose VEV Xt(o) '■= {0\X t (q)\0} is also known as the topological susceptibility of QCD 2 , 
are presented here. This correlator has been connected to the mass of the r\ -meson through 
the Witten-Veneziano formula [8-11]: 



« 2 

32ivr 2 



Xt(q) 



(leading order) (1-10) 



11, 



where F n ~ 94 MeV is the pion decay constant. An explicit sum rule calculation with an 
OPE at one-loop level using a Borel transformation has been done in [12]. In this work the 
value m 1 ~ 1 GeV is correctly estimated. 

The correlator defined in (1.1) with renormalized operators is finite, i.e. all its matrix ele- 
ments are finite, except for possible contact terms. These arise from the point where x = 



Effectively this expansion separates the high energy physics, which is contained in the Wilson coef- 
ficients, from the low energy physics which is taken into account by the VEVs of the local operators, the 
so-called condensates [4]. These cannot be calculated in perturbation theory, but need to be derived from 
low energy theorems or be calculated on the lattice. 

2 For a discussion of topological effects in QCD and the significance of the operator Oi and the correlator 
(1.9) in that respect see e.g. [6, 7]. 



and manifest themselves as divergences oc 5{x) and derivatives of 5(x) or in momentum 
space terms polynomial in q. These local terms do not contribute to sum rules and can 
and should be subtracted with proper counterterms. 

The leading term on the rhs of (1.3) is the coefficient in front of the unit operator 1 which 
is just the perturbative VEV of the correlator (1.1): 

(Q 2 ) 2 C (q) = (0\X t (q)\0)\ pen . (1.11) 

The coefficient Co is known for the scalar case (1.5) at four-loop level [13] and for the 
pseudoscalar case (1.6) [1] and the energy-momentum tensor correlator [2] at three-loop 
level. The next important contribution in the OPE is the coefficient of the dimension 
four operator [Oi] (1.5). 3 The coefficient C\ has been calculated at two-loop level for the 
scalar 4 and tensor cases [2]. Here we present the coefhcent C\ for the pseudoscalar case at 
three- loop level which so far has only been known to one-loop accuracy [16, 17]. 

All physical matrix elements of [0\] = ZqO^ are finite and so is the renormalized coefficient 
d: 5 

d = -^-Cf. (1.12) 



Z G 



The renormalization constant 



aB _l + ._ln^-(l-£LJ) (1.13) 

has been derived in a simple way in [18] (see also an earlier work [19]). Here Z ag is the 
renormalization constant 6 for a s and the /3-function is defined as 

d , ^„ /cO m 

i>0 



P{a 3 )=^ 2 — 1 \na a = -Y J fii[ — ) • (1-14) 



The outline of this paper is as follows. In the next section the renormalization properties 
of 0\ will be discussed. In section 3 the details of the calculation will be described (section 
3.1) and the results for the OPE of (1.9) will be presented (section 3.2). After that Renor- 
malization Group invariant (RGI) operators and Wilson coefficients will be constructed 
(section 3.3) followed by a numerical evaluation of the main results (section 3.4). Finally, 
some conclusions and acknowledgements will be given. 



3 In the case of massive fermion flavours / we would also have contributions proportional to the dimension 
two operator O* = m 2 f 1 and the dimension four operator 2 = mfipfipf. In the case of temperature T^O 
Lorentz variant operators like Tq ~ e + p with the energy density e and the pressure p have to be considered 
as well. At T — 0, however, only Lorentz and gauge invariant scalar operators contribute to the the VEV 
in (1.2) which is the quantity that we are ultimately interested in. For a discussion of the correlator X t (q) 
at finite temperature up to 0(a B ) see [14]. 

4 The one-loop result for the scalar case was first derived in [15]. 

5 In the massless case 0\ only mixes with unphysical operators whose matrix elements with physical 
external states vanish. The renormalization of 0\ including these unphysical contributions as well as the 
mixing with 2 in the massive case can be found in [18]. 

6 Often in the literature Z aa is used instead of Zg and a^C" G M „ instead of 0%. This renormalization 
is only valid up to first order in a s as the renormalization constants Zq and Z as coincide to this accuracy. 
In higher orders, however, Zg and Z as differ. 



2 Renormalization of 0\ and its correlator 

The operator 0\ forms a closed set under renormalization with the pseudoscalar fermionic 
operator 

d,4 ■■= £ ™^£* /7w 7„ 2 7„3*- (2-1) 

/ 

which can be written as 

-Vs'^E^/W (2.2) 

/ 

in the Larin scheme for 75 [20]. 

The e-tensors appearing in (1.6) and (2.1) are then drawn out of the R-operation per- 
formed in dimensional regularization. In the correlators which have to be calculated there 
are always two e-tensors involved which can be contracted and expressed through metric 
tensors: 

P MlM2M3M4p _ _ n [Ml n A»2 „M3 „M4] o o\ 

t t -V 1 V2lS3V4 ~ il V\ii V2i> vzil 1/4 l Z -°7 

where [. . .] means complete antisymmetrization. These operators are renormalized like [20] 



[d,4] = ZiZUsW? 11 = Z s jd,J^, (2.4) 

[d l ] = z GG df + z GJ d^i» (2.5) 

where Z S MS is an MS renormalization constant, Z| a finite renormalization constant fixed 
by the requirement that the one-loop character of the axial anomaly relation 

[d,4\ = ^n ! T F [O l ] + CT (2.6) 

is valid in dimensional regularization. 7 CT stands for contact terms of d^J^ with fermion 
fields. In the gluon sector these can be neglected. Z G q is an MS renormalization constant 
again and Zqj starts at 0{a s ). In [20] Z S MS and Z| are given up to 0{a e ) and 0(a 2 s ) re- 
spectively. Furthermore it is shown that Z GG = Z a (Z a being the renormalization constant 
for a 3 ). The constant Zqj is only given at one-loop level in the literature [1, 20] but for the 
Wilson coefficient C\ at three-loop level it is needed to two-loop accuracy. In section 3.3 
we will also need the corresponding three-loop anomalous dimension. The simplest way to 
determine Zqj is by constructing the matrix elements of 0\ and d^ J 5 with two external 
fermions (see Fig. 1) using a projector 

p(g):= ^y* V 3 7 M4 w ^ 4 (27) 

on the external fermion line. From this we get 



7 In Pauli-Villars regularization for example this relation is automatically fulfilled. In d 7^ 4 dimensions, 
however, the operators d^J^ and Ox become linearly independent. 
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Figure 1. Diagrams for the calculation of Zqj 
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(2.8) 



y y 

An interesting additional application of this result is to check the connection between the 
anomalous dimensions of the operator set {Oi,d u J^}. In [20] the following relations have 
been motivated: 



/%.) 



with 



7ij 





7GG 


( 




r 1 






1GJ 


"( 


4^', 


) 


7j 


/' : 


d 7 " 


K* 




z 


-I 



Z G g Z GJ 

z s r 



(2.9) 
(2.10) 

(2.11) 



The first relation (2.9) has been explicitly checked to three-loop level in [20] the second one 
(2.10) only to one-loop accuracy. Now we can check this equation with "fGJ a t two-loop 
level and 7j at three-loop level and it turns out to hold there as well. Using (2.8) and 
the renormalization constants Zj and Z a [20, 21] the following anomalous dimension is 



derived: 8 

7GJ = -12C, (J) + (^) 2 {-^C^C, + 36C| + fc F n f T F 

+ ^y[-^fc 2 A C F + A6lC A C F + 576C A C F n f T F C3-^fc A C F n f T F (2.12) 

-126C F - 576C 2 F n f T F ( 3 + 428C 2 F n f T F + ^C F n 2 T 2 

6 

Now we can write the correlator Xt(q) as 

i [d^xe^Tild^ld^O)} 

^^{zjgof (x)6f (o) + 2Z G(5 z GJ of (*)0 m j?"(o) + ;&,M? "(*)a,jf *(o)} . 

(2.13) 

In [2] it has been discovered that there are contact terms at two-loop level in the coefficient 
C\ for the correlator of 0\ . The coefficient C§ also has contact terms for the correlator of 
two operators 0\ or two operators T^ u . For the operator 0\ we can make an important 
restriction on possible contact terms due to the fact that it can be exactly expressed as the 
divergence of the Chern-Simons current: 

6\ = d„K^ (2.14) 

with 

K* = e» vpa \AG a v d p G a a + -gJ ahc GlG b p Gl J . (2.15) 

From this follows for (2.13) 

.yVze^TtidiKxNdiKo)} 

=q^iJd 4 xe^T{z 2 Gd K B ^x)K B ^0) + 2Z Gd Z G jK B ^x)J^(0) + Z 2 G jJ^(x)J^(0)} 

^q,q» {<WV) + ^CfV) + • • •} for q 2 -». -oo (OPE) 

(2.16) 

with dimensionless coefficients Cq U {q 2 ) and Cf (g 2 ). Because of the non-local factor -% 
the coefficient C^ 1, (q 2 ) cannot contain any contact terms. This makes the Wilson coefficent 
G\{q 2 ) = -^-r-Ci" (q 2 ) for the correlator (2.13) finite and unambiguous due to the absence 
of contact terms. 
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Figure 2. Diagrams for the calculation of the coefficient Cq(Q 2 



3 Calculation and results 



3.1 Details of the calculation 



The leading coefficient Cq is just the perturbative VEV of the correlator eq. (2.13) 



(QTO)(g) = (ON /d^e^TflOiKaOPiKOtflO) 



port 



(3.1) 



which has been computed up to order a 2 (three loops). In Figure (2) some sample Feynman 
diagrams contributing to this calculation are shown. The operators Of and d^J 5 ^ play the 
roles of external currents. The Feynman diagrams have been produced with the program 
QGRAF [22]. As all diagrams in this problem are propagator-like the relevant integrals 
can be computed with the FORM package MINCER [23-25]. For the colour part of the 
diagrams the FORM package COLOR [26] has been used. 

In order to compute the coefficient C\{Q 2 ) the method of projectors [27, 28] has been 
applied, which allows to express coefficient functions for any OPE of two operators in 
terms of massless propagator type diagrams only. The method is based the fact that in 
dimensional regularization every massless tadpole-like Feynman integral is set to zero. 
We apply a projector to both sides of (1.3) which sets every operator on the rhs to zero 



8 7gg an d 7j can f> e found in [20, 21] at three-loop level. All renormalization constants and anomalous di- 
mensions are available at http://www-ttp.particle.uni-karlsruhe.de/Progdata/ttpl3/ttpl_3-003/ 



except for 0{ 



B. 



r{x t ( q) }=j:m 4 -^^c^mi>{o?} 



(3.2) 



with P{0^} = 1 and P-fO^} = 0. This is done in the same way as described in [2] 
leading to 



,2 r,(6?,6f) (n2 . 



with 



(of ,a„j 6 -) (g2) + zhc ^ >dv jB* )(n ^ 



'(Q 2 



(3.3) 



cT°\q 2 ) 



5 ab g^ 2 1 d d 
~ g (D - l)Ddh"dki 




9b fi 2 



, (3.4) 



ki=0 



where the blue circle represents the the sum of all (bare) Feynman diagrams which become 
1PI after formal gluing of the two external lines representing the operators on the lhs of 
the OPE. 
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IQB qB 

Figure 3. Diagrams for the calculation of C\ g ' 1 
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d v J, Bv 




+ ... 



(O b d J B ") 

Figure 4. Diagrams for the calculation of C\ ^ ' " 5 . 



Table (1) shows the number of diagrams generated for the different contributions to Cq 
and C\ All results are given in the MS scheme with a s = ^r, a s = |i and the abbreviation 
l^ q = In ( %x ) where [i is the MS renormalization scale. They can be retrieved from 
http://www-ttp.particle.uni-karlsruhe.de/Progdata/ttpl3/ttpl3-003/ 



d u J, BV 
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duJt 




+ • 



(d J B " d J 1 

Figure 5. Diagrams for the calculation of C\ g 5 " 5 



Correlator 


loop 


1 loop 


2 loop 


3 loop 


(0|Of(x)Of(0)|0) pert 





1 


12 


215 


(o\6f(x)d v j^(o)\o) peIt 








1 




{0\d v J^(x)d 1/ J^(0)\0) peTt 





1 






Px(6f(a;)t^(0)) 


2 


75 


2567 


94964 


Pi(O?(x)d v J?"(0)) 





8 


345 




Vx{d v j£ v (x)d„j£ v W) 





8 







Table 1 . Number of diagrams needed for Cq and C\ 



The gauge group factors are defined in the usual way: C F and C A are the quadratic Casimir 
operators of the quark and the adjoint representation of the corresponding Lie algebra, d R 
is the dimension of the quark representation, n g is the number of gluons (dimension of the 
adjoint representation), T F is defined so that T F 5 = Tr (T a T J is the trace of two group 
generators of the quark representation. 9 For QCD (colour gauge group SU(3)) we have 
C F = 4/3 , C A = 3 , T F = 1/2 and d R = 3. By n f we denote the number of active quark 
flavours. 



3.2 Results 

As we have seen from (2.13) contact terms in Co are possible and it turns out that they 
appear starting from one loop. Because of these contact terms an unambiguous result 
for Cq can only be given up to local (that is q-independent) contributions. To avoid the 



9 For an SU(JV) gauge group these are d R = N, C A = 2T F N and CV =T P (N - i) 



ambiguity the Q 2 -derivative is presented: 



2 d r _ n 9 



-1 + aA ~t^Ca + -zrifTp. j + a s l m ( — — C A + -n f T F j 

2 / 51959 2 107 3793 251 2 2 

+a C , H nfl F C F H n,C A l F re.i „ 

3 V 864 A 12 ' 108 s 54 / F 

55 1 \ 

+ yC 3 Cl - Kzn s T F C F + -Can^xT^J (3.5) 

+a% q {~Cfl A + 2n / T F C F + fn f C A T F - 7 -n^ 

This result has been derived before [1] which serves as a nice check for the setup. As 
discussed above the coefficient C\ is unambiguous and is therefore given in full: 

d =64 1 1 + a s f "og-C^ - n n / T f + ToU^ ~ a"U n / r H 

2 /25945^ 2 11 m ^ 4355 ^ m 25 2m2 1727, <, 

+a 2 C 2 n f T F C F n f C A T F + —n 2 f Tl + L C 2 

s V 1296 A 2 f 648 s 81 ' F 216 "* A 



i nTC -Mi n CT +™l n^' + ^C 2 -^ 



1 33 3 

+5<™"? T * - y&cj + 3C 3 .!.,r F C F - -C 3 n,C„T F 



/ 19360399 461 _ 614501 

s V 186624 A 144 7 F 10368 / 



1857805 - 28981 2 2 126415 2 2 125 3 3 
31104 / A 2592 / F 15552 ' F 729 ' F 

+»=»i„cl + ?W,cj - l ™i„ q „,c A T r c r - ^i„,n,c°T F 

10368 m A 32 M9 / F 64 m<7 / a f f 1?28 M9 / A 

105 2-2^ 6661 , 2^ rj.2 25, 3^3 9779 2 3 
+ T6~ l M n f F F + "goT ^ n f A F ~ gi l M n f T F + "864" v?^ a 

_^/ 2 „ r t r - 2795 / 2 47 r 2 r 4- — / 2 T7 2 r 2 r 4- — / 2 47V t 2 

96 V 9 n /°AJF^F 288 ^ q n f u A i F + ^i,,"/^^ + 24"' f -* F 

Q Q Q lOOl O Q A u J_ _ Q O 1 J. Q Q ^^ Q 

~^fK q n f T F + yj^ l M C A - Y^ l M n S C A T F + 7^Ki n f C A T F 

1 55 15 

-^ft n3 f T F + yCsCi - 15C 5 %T F C| + -( 5 n f C A T F C F + 5( 5 n f C 2 A T F 

6893 3 145 ^1291 349 2 

~l44~^ 3 ^ + J2^ 3Uf F F -^Csn f C A T F C F - —(3n f C A T F 

— 7r(3 n fT F C F + -^-(3n f C A T F — ■^rC^m^A + -rC3^ q n fC A T F C F 

2 J 36 32 4 



—3(3l M n f T F C F + -Qil M n f C A T' F 

The cancellation of all divergences is a strong check for this result. Another important 
check is the independence of the gauge parameter £ as all calculations have been done for 
an arbitrary R^ gauge. The leading term of (3.6) is in agreement with [12] and the part 



10 



oc a s / Mq has been derived in [17] if we set the colour factors to their QCD values 
QCD colour factors we get 

C] =64 < 1 + a, ( n f -\ l,.„ l,.„n« \ 

\ V 12 18 / 4 m 6 M * V 



10 



For 



+a 



2 / 25945 5939 

^ 144 ~ 432 



25 

n. 

324 f 



n f + T^-rn, H ^—l 



1727 



62 5 2 121 2 

-yU n / + 54 W 1 / + -jjU 

11 2 1,2 2 297 

12 M « y 36 m f 8 



rC3 n 



+a* 



/ 19360399 7972411 



V 6912 



20736 



-n, + 



611093 2 125 o 

n, n f 

62208 ' 5832 7 



594247 

384 

9779 o 

1 

32 M9 



25 



264113 9181 2 _ 

"TT52~ M% + 1152 mU i ~ 648 "'"' 



9485 



192 



■*>/ + 



649 



;2 2 



216 



Z 2 n 3 



+ 



1331 



64 



121 



11 

48 M9 



t n n s + ^/Ln 2 



1 



1485 A 145 . 
— C 5 + -^5% 



—J 3 r> 3 

216"" 

20679 . 46333 , 

C3 + OOA C3% 



16 '" 864 
17. 2 9801 , 33 A , 1 , 2 

+ 48^ 3n -f ~ ~32~^ 3 M + T^ n/ + g^^ n / 

A nice consistency check for these results is to perform an OPE of the correlator 

i J ^xe^T{[d^]{x)[d^m} = {Q 2 fC{ J + d JJ [Oi] + . . . 
and then see that (2.6) is fulfilled (except for possible contact terms): 



d_ 



Q 2 ^Co J =(^n f T F r Q ^C , 



C{ 



jj 



• Ik 



Indeed we find 



Q 2 ^C J n J - Ha 



and 



iJJ 



,2 m 2m2 



C( J =4ainiT£ 1 + a s [ —C A - -n s T F + 



dQ 



/157 



2 ^0 



7T- 



n f T F ) 2 d. 



-^n 2 T 2 F 
16 > F 



V 36 



9 



11 
12 



■p.q^A 



1 



;l M nfT F 



(3.7) 



(3.8) 

(3.9) 
(3.10) 

(3.11) 
(3.12) 



satisfying (3.9) and (3.10) up to the calculated accuracy of 0(a 2 ) and 0(a 3 ) respectively. 



10 In [17], however, the leading term differs from this result and the one derived in [12] by a minus sign 
and the non-logarithmic term of 0(a s ) is also missing there. 
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3.3 RGI operators and Wilson coefficients 

Note that the coefficients (3.5) and (3.6) are not Renormalization Group invariant (RGI). 
In this section we take RGI versions of all operators and construct RGI Wilson coefficients. 
For an operator that is renormalized multiplicatively like d^J^ in (2.5) constructing a finite 
and RGI operator is straightforward (see e.g. [29]). Because of 

^MI = -#(«.(/*))[«] ( 3 - 13 ) 

we can define 

/ ^B) da M (3 - i4) 



=:E 2 (a s ) 

which fulfills /i 2 ^-[<9^ J^] RGl = 0. A remarkable feature of the operator (3.14) is its renor- 
malization scheme independence (see also [30]). If we start with a different renormalized 
operator 

[d,4\' := Z{a s )[d,4\ (3.15) 

we get 

-fjifls) = YAa s ) + fi 2 -^ HZ(a s )) (3.16) 

which leads to 

„i i v E 2 (a s ) . . 

E 2 {a s ) = (3.17) 

Z(a s ) 

and therefore to the same RGI operator 

[«V 5 T GI = E' 2 (a s )[d,4]' = E 2 {a s )[d,4). (3.18) 

If we apply the same procedure to the non-diagonal operator 0\ we get an RG variant 
operator 

a s (p) 

f 7gg( q ) 



i$U da * [c 



=:Ei{a s ) 
where E\(a s ) = a s because of (2.9). Taking the derivative wrt the renormalization scale 
we find 

/^[Oi] RGV = E 1 (a s h GJ (a s )[d^} = §^}fGj(a s M4f G1 (3-20) 

which leads to the definition of the RGI operator 

a 3 {p) 

\n iRGi \n iRGV f Ei(a) , . da ^-.rgi 

J E 2 (a) a (3(a) ^ 

=:a B Z(a s ) 

=a s {z GG (a s )6f + (z GJ (a s ) - E 2 (a s )Z(a s )Z s j(a s )) d^J^) 
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fulfilling /i 2 ^pj[Oi] RGI = 0. In similar way as for (3.14) it can be shown that (3.21) is 
invariant under transformations [Ox] — > [Ox]' = Zx(o s )[Oi]. Even if we allow for redefinitons 
of the kind [Ox] -> [Ox]' = Z 1 {a 8 )[dx] + Z 2 (a s )[d^ Jg] the RGI operator derived with this 
method is the same: 



[d x ] RGV/ = [0!] RGV + |i^4l4^ I^T GI 



^^-r^[o 



=S iRGV' 



dfJ, 2 



E 2 (a s )Z 1 (a s ) 

Ex(a s ) , . 2 d /£ , i(a s )Z 2 (a s ) 
7Gj(a s ) + At Vo 



-M^ 



[O: 



]RGI' = [(5i] RGV 



a 3 (ji) 



dfi 2 \E 2 (a s )Z 1 (a s ) 

Ex(a) da 

^2 (a) 



Ml 



RGI 



(3.22) 
(3.23) 



+ 



Ex(a s )Z 2 {a s 



a(3{a) I £ , 2 (a s )Zi(a s ) 



M 



RGI 



=a s Z'(a 3 ) 

[dxf GV -a s Z(a s )[d,4}* G1 = [dx} KG1 . 



-RGI 



(Q 2 ) 



2\2 



(0\Xr\q)\0) 



pert 



(3.24) 



(3.25) 



(3.26) 



can now be calculated from the same three bare correlators as Co and the result for its 
Q 2 -derivative is 



The leading RGI Wilson coefficient 
in an OPE of the RGI correlator 



X^ Gl (q) :=i d A xe iqx T{[dx] RG \x)[6x] RG1 (0)} 



Q 



d (jRGl _ a s n g 



dQ 



2 ^0 



TT Z 



+ 



97 7 11 2 



{\1C A - 4n f T F 



■18n f T F C F + a 2 s 



51959 2 107 



864 
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3793 251 22 1135 2 

+ 108 n/ 5X n/ 48" 4 ~*~ M?n/ 

46 7 2 2 121 2 2 11 2 

+-^l m n f C A T F — -l M n f T F tjT^C^ + "^"V^/Ca^V 



(3.27) 



+ 



(HC^-^Tjr) V 2 



/291 on 

-n f C A T F C F - 42n)TlC F 



2r7-i2 / 



-\-66l^ q nfLy A l F U F — IZl^qTbfl F C/jr 



+ 



297 



- 4n,T F )2 T ^^ 4 — ^/OVCV 



475 
T 



-108n 2 T|C| - 37n 2 (7 7l Tj<7 F + 4^1 CV 
An explicit calculation confirms that indeed H 2 -^z \Q 2 ~d(ji O rgi \ = 0- 
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275 
"96 



(3.28) 



As explained in [2] a finite and RGI version of 0\ can be denned as 

Of GI := $(a.) [Oi], ${a.) := ^T^ = a s ( 1 + £ fa 
The RGI Wilson coefficient 

P\P>s) *■ 

+(2Z G(5 Z GJ - 2£ 2 Z Gc5 Z5Z)C7J f AJ "^(Q 2 ) (3-29) 

+ (^ - 2£ 2 Z GJ Z5Z + (£ 2 ZjZ) 2 )c{ 9 /* ^^(Q 2 )} 

which satisfies 

Cf Gl [O l \ KGl = C l [O l \ (3.30) 

in the OPE of (3.26). The result is 

C? GI =64a s |l + a s (^C A - ^n f T F + ^l M C A - \l M n f T F \ (3.31) 

+ (nc A -\n f T F ) (-T°2 " l^T^, + ^T,) 

2 ^ 25945 ^ 2 11 4355 25 2 2 1727 2 

+ ° s vT^e -6 ^ " T n ^ c ' " "61^ %6aTf + 8i n ^ + "216" UCjl 

limThfJ. pOp ~" £ ,,nTh fU 4 -i p — T — ~ I ,, n Tl n-L m ~~ r" fc O * £ Tl f\^/ A ± J? 

r\ HQ J -r r c\f-7 A"2 J A i* ' qj-t (*<i' v f F ' -i A A fj,q A -| n /^g J A r 

a * f 2669 C 3 » nCTC »„ c2T 

+ (11^-4^^) v " 72 Ca ~ 12 n f C *l F C F + 36 n f C A l F 

+— n 2 r 2 r - —n 2 r t 2 - —i r 3 - — / n r r r + — / n r 2 r 

t o n f i F ^ F n f u A i F 4 m 9 u a x V8 n / u A J iJ-^F -r -,„ M g 7 A F 

+5/ n 2 T 2 C - -Z n 2 r T 2 ^i 4- a ' /10619 4 

+ ^n C T C 2 + 1 -^n C 2 T C + *°H„ C 3 T + ™n 2 T 2 C 2 
t n f u A i F Ls F -t- n f ^ A i F ^ F -f %W i f + 9 n f 1 F u F 

7r~ nfC A l F C F — rij:G A l F — —njl F C F + ~To n f A f ) 

/ 19360399 461 _ 614501 

s V 186624 A 144 7 F 10368 ' 

1857805 o m 28981 ,„- 126415 2 ^ m2 125 , m » 

n f C 2 T F + n 2 T 2 C F + n 2 C A T 2 n 3 T 3 

31104 J A 2592 ; F 15552 f F 729 ' F 

594247 35 _ IM^ _ JW^*,. 

10368 M9 A 32 ^ 9 J F 64 M9 ' 1728 " q J A F 

105 2™2^ 6661 7 2^ ^2 25 33 9779 2 3 

+ Te" ^ n f F F + ~8Q4~ >" qn f A F ~ ~£l l M n f T F + "864" V«^ 
275 / 2 n rTr - 2795 / 2 „ r 2 r 4- — / 2 r7 2 T 2 r 4 — / 2 « 2 r T 2 

l ^q n f L/ A^F^F 2gO L M n f i ~'A 1 F "I" 2A M f F^F "t" 2^ l m U i^ A->- F 
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-SCaUn^C, + ~ 2 C^y s C A T^ + ^ -^ 



O O Q ^^ Q lOOl Q O X^Xo Q ^^ J. J. o r> ^^ r\ 

-^Ki n 5 T F + yj^ l M C A ~ Y^ l m n f C A T F + -^ l ^q n f C AT F 

1 55 15 

Yft n3 f T F + y CsCl - 15C 5 %r F C| + yCsnyC^T^C^ + SCsn.Cir, 

6893 3 145 2 1291 349 2 

-^4jC3C A + — ( 3 n f T F C F + ^-C3^/C A T F C F - — ( 3 n f C A T F 

-^-C^n f T F C F + -rT-Cz n fC A T F — ^r(3l M C A + —C, 3 l^ q n s C A T F C F 
I ' 3d 62 4 

,2^ , 3 A7 „ 2 ^^ , o? f 441065 4 

2592 A 

109529 2 1415 3 165 2 2 2 15835 2 2 

432 n/ "' n - n/ A F "' n - n/ F F + 216 n/ 

7825 2 2 2 125 33 125 3 3 29359 4 

~~2i6~ n/ A F " ^r n / T ^ + ^r n / 6 ^ - i^* 6 * 

7717 2 5281 3 45 2 2 2 925 2 2 

j^~lM n fC A T F C F H ——l m n s C A T F + —l m nfT F C F + -—l M n f C A T F C F 

W~^q n f^' A ^F W'^q n f^F^F + ~< M9 ^ /C/ A Tp 988 M9 

-^>^ F C F + ^>^T F + f «C A r|C F - 4*>^2 T| 
-| W3?C7, + jj^n?^ + ^C\ + ^C 3 %CiT F C F 

-f Cs^T, - 45fcn*Z*C? + ^ C 3^C A TjC F - yCs^T, 2 

of / 1667183 ^ 5 29359 2 ^2 227807 3 

+ (llC A -4 % r P )2 v ~^036^ + -96~ n/C ^ J ' C ' + "lW"' ^^ 

1525313 4 727 2 2 2 33923 2 2 2 129511 2 3 2 

"! 5184 Hf A + "3" n f C A T F C F ^^-n f C A T p C F ^j— n f C A T F 

215 3^^ , HZ 3r7 T 3 r , 10949 „3 r 2 T 3 , A ?? 4 r 4 r - — « 4 r T 4 

n f i F i^ F -\- n / u A J F o F -t- n f u A i F -+- ?i f i F ^ F \r<i s 

116809 f C 5 + 2057 nC 2 Tc2 «! nC 3 TC 120667 4t 

3456 ™ * + 09 'cr/^-t-'^FT _„„ ^ q n f ^ A i F ^ F ^- L M n f L, A ± F 

I 143 / „2 r T 2 r 2 _ 897 / 2 r 2 T 2 r _ 12989, „2 r 3 T 2 _ f^, ^t^^ 
"I" . t M9 'l / '-'A-t F ( -' F j_ l m a j Xy A 1 F^F 200 i w' i f^A 1 F 9 t M9' i / J F Ly F 

197 / „3 r T 3 r , 346 > „3 r 2 T 3 , 1, „4 T 4 r _ J_^ 7 4 r T 4 

t- 1 m n s k_j A ± F <^ F -t- i m n^ A i F -t- L^ q n f i F i^ F — i m n } ^ A ± F 

I "* f_Z^3_ n C 2 T C 3 22121 3 2 31207 4 T c 

+ (llC A -4n/T F )3 I 16 n/GAi '°* ,+ 32 n /^^ + 32 rt f C A l F C F 

ZU I t/ Q Q q J.OOOO o O ^^ Q o ^i t? O ^ I O Q O Q n Q 

- rij:G A l F C F H - n f G A l F C F — n^U A l F G F — 45n^i F C F 

-^n^C A T|C F + 1443n3c*lT|C F + 178ra 4 T F C F - 384n 4 C A T 4 C F + ^1^T|C F 
-2178C 3 ^Clr|C| + 2178C 3 ^C3T|C F + 1584C 3 n3C A T|C F - 1584C 3 ^CJT|C F 
-288C 3 ^T 4 C F + 288C 3 n 4 C A T 4 C F )} . 

Again an explicit calculation confirms that indeed fi 2 -j-^Cf- G1 = 0. This result can now 
be used to obtain the logarithmic pieces of Q 2 ^n^Co lGI an d Cf- G1 at four-loop level. If a 
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generic RGI quantitiy has the structure 

Q RGI = a s (ix)A l + a s (fi) 2 (A 2 + l m B 2 ) + a s (^f(A 3 + l, q B 3 + l 2 M C 3 ) 

+ a s ^)\A A + l m B A + l 2 m C A + ll q D A ) (3.32) 



f(A 5 + l m B 5 + ll q C 5 + ll q D b + l^E 5 ) + 0(o«) 



with scale independent coefficients (j4j, Bj, . . . ) the requirement ^ 2 xt2"Q RGI = leads to 
the conditions 

B 2 = Ai/9b, 

C 3 = B 2 /3 , J B 3 = A 1 /3 1 + 2A 2 /3 , 

D 4 = C 3 p , C A = X - (3B 3 f3 + 25 2 /3i), 5 4 = Ai/9 2 + 2A 2 /3x + 3A 3 /3 

which in the cases of Q 2 ~^iCq and C 4 can be used as checks and 

E 5 = D A (3 , 

As = 5 (4C 4 A) + 3C 3 ft) , 

C 5 = l - (2B 2 p 2 + 3B 3 0i + 45 4 /? ) , 
B 5 = A ± /3 3 + 2A 2 /3 2 + 3A 3 ft + 4,4 4 /3 . 

Using the four-loop /3-function 11 of QCD [37, 38] the following four-loop contributions (for 
QCD colour factors) are derived: 



i d rgi, 4ioo P «X ![!!/_ HlS 121n / 1331 \ 3 
ciQ 2 ° vr 2 1 I 54 12 + 8 16 ) M 



/7n3 1783?i 2 21647n / 19569 \ 2 
+ ^"36" 144 + 96 16/ " 9 

1 /251n« lOnfe 147169n* (3.33) 

+ (33-2 % ) 2 l^l" + ~ 432" " 330n /<3 



108663n^ „ 48109321< 299475n,0 

+ i + 10890n 2 C3 f - — 

8 /S 192 2 



138470387™ ,■ 5929605C 3 450379545 \ , 

1 (,,„ + const. 

64 8 64 J m 



L The one-loop, two-loop and three-loop results are known from [21, 31-36]. 
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ci 



RGI,41oop = 64a 5 



7? 



/ 



1296 216 



lln 3 t I2ln 2 f 



1331n f 14641 \ 






96 



+ 



256 



V 



+ 



'hn) 1595^ 4355n 2 293975% 424105 \ 3 
972 ~ 2592 + 216 1152 h 384 J " 



(ihn) n 3 f ( 3 6937^ 77n^3 1812625™ 2 



1944 



24 



2304 



16 



13824 



617171^3 95413377., 10781 1( 3 12658057 \ 2 
+ 32 512 64 + 1536 J Mg 



+ 



1 



125n^ 17n^ 3 457613n^ 4237n^3 



(33 - 2n f ) 2 \ 2187 18 15552 

580ra^5 1320687777^ 38583n^3 



54 



9 



+ 



5184 



+ 



+ 2695n;C5 



+ 



90573423577^ 1172479n 2 -C 3 56265n^5 6551159345n 2 

10368 4 2 + 4608 

11374907577^3 45919577^5 16816549087n / 

32 4 1536 

486694791C3 17788815C 5 48864828943 \ 



+ 



+ 



L„ + const. 



32 8 1536 

For completeness we also give the RGI Wilson coefficients for the correlator 



/ dW^r^j^^Ml^o)} = (Q 2 ) 2 c JJ ' RGI 



The results read 



+ C( J ' RGl [Oi 



RGI 



+ 



cy 



and 



C{ 



JJ.RGI 



dQ 2 ^ 
7 157 



JJ,RGI n g 



TT Z 



-^n 2 T 2 
16 s F 



, s n 2 T 2 ^ l + a s i — ^ 



C A 



5 111 

-n s T F + — l m C A — -l^rifTp 



+- 



(Is 
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C 2 A - 15n f T F C F + 5n f C A T F j j 



(11C A - 4n f T F ) V 2 
The four-loop extension of these results with QCD colour factors are given by 

n 2 .(-33 + 2n / ) 



q2 d qJJ,RGI, 41oop _a s n g 



dQ 



2 w 



7T' 



384 



+ const. 



and 



C JJ,RGI, 41oop = 4a 3 f , 1 n 2 (1416 g _ 1533n + 20 2) 

L 864 



+/ 



121nj 
64 



Unj n) 



18 



144 



+ const. 



(3.34) 



(3.35) 



(3.36) 



(3.37) 



(3.38) 



(3.39) 
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3.4 Numerics 

We now consider the two cases n f = (pure gluodynamics) and n f = 3 which are most 
important for applications. Furthermore we set Q 2 = fi 2 , i.e. l m = 0. The numerical 
results for C\ and Cf 1 are then 

d(Q 2 = //, n s = 0) = 64{1 + 13.0833a s + 135.547a 2 + 1439.88a 3 }, (3.40) 

d(Q 2 = n 2 , n s = 3) = 64{1 + 12.25a s + 94.0971a 2 + 646.69a 3 }, (3.41) 

Cf Gl {Q 2 = /i 2 , n f = 0) = 64a s {l + 10.7652a s + 102.475a 2 + 1089.78a;?}, (3.42) 

Cf Gl {Q 2 = /i 2 ,% = 3) = 64a s {l + 9.13889a s + 55.9532a 2 + 361.615a 3 }. (3.43) 

In order to estimate the numerical significance of the higher order corrections we evaluate 
C x at n = M z , n = 3.5 GeV and fj, = 2 GeV with 

a^f =5) (M z ) ~ 0.118 , a^' =3) (3.5GeV) w 0.31 and a( n / =3 )(2GeV) w 0.47 [39] (3.44) 

for the cases n f = 5 and n f = 3 respectively. 

Ci(Q 2 = i? = M z ,n f = 5) = 64(0.0116 + 0.0949 + 0.4393+1), (3.45) 

3 loop 2 loop 1 loop 

d(Q 2 = i? = (3.5 GeV) 2 , n s = 3) = 64 (0.6213 + 0.9162 + 1.2088 +1), (3.46) 

3 loop 2 loop 1 loop 

diQ 2 = n 2 = (2 GeV) 2 ,^ = 3) = 64(2.1654 + 2.1061 + 1.8327+1). (3.47) 

3 loop 2 loop 1 loop 

At the scale /U 2 = Mz the two and three-loop contributions are about 9% and 1% wrt 
tree-level, whereas at a scale /z 2 = (2 GeV) 2 these contributions become so large that 
perturbation theory stops to work (as is expected). From this evaluation we can assume 
that in the case of Q 2 = \j? the Wilson coefficient to this accuracy in perturbation theory 
is a valid approximation down to a scale of about fi 2 = (3.5 GeV) 2 . 

4 Discussion and Conclusions 

I have presented higher order corrections for the coefficient function C\ of the OPE of the 
correlator Xt of two pseudoscalar gluonium operators. This result extends the previously 
known accuracy by two loops. It is also worth of notice that no contact terms can appear in 
this coefficient due to the relation between the operator 0\ and the Chern-Simons current, 
a fact that has been explicitly checked and verified up to O(a^) by this calculation. The 
OPE of the correlator of two operators d^ J§ which mixes with 0\ under renormalization 
has been performed as well and the corresponding coefficients Cq and C{ have been given 
at three-loop level. In addition the construction of RGI operators and Wilson coefficients 
has been discussed, the coefficients C^ Gl , Cf 1 , C ' and C 1 ' have been presented 
and their logarithmic part has been derived at four-loop level from the principle of scale 
invariance. Finally, a numerical evaluation shows the validity of the OPE for the important 
case n f = 3 and the choice Q 2 = fi 2 down to a scale fi 2 = (3.5 GeV) 2 . 
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